High-precision determination of transition amplitudes of principal transitions in Cs 

from van der Waals coefficient C 6 



Andrei Derevianko 

Department of Physics, University of Nevada, Reno, Nevada 89557 
(Dated: February 2, 2008) 

A method for determination of atomic dipole matrix elements of principal transitions from the 
value of dispersion coefficient C§ of molecular potentials correlating to two ground-state atoms is 
proposed. The method is illustrated on atomic Cs using Cs deduced from high-resolution Fes- 
hbach spectroscopy. The following reduced matrix elements are determined (6SV/2II-DH6P1/2) = 
4.5028(60) |e|a and (6S1/2II.DII6P3/2) = 6.3373(84) |e|o (o = 0.529177 x 10~ 8 cm.) These matrix 
elements are consistent with the results of the most accurate direct lifetime measurements and have 
a similar uncertainty. It is argued that the uncertainty can be considerably reduced as the coefficient 
Cs is constrained further. 
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The leading long-range interaction of two atoms in 
their respective spherically-symmetric ground states is 
described in terms of van der Waals coefficient C® [? 
] . Studies of magnetic- field induced Feshbach resonances 
and photoassociation in ultracold atomic samples allow 
to constrain this coefficient. Several highly-accurate de- 
terminations of Cg for alkali-metal dimers were recently 
reported [?????]. Here we propose a method 
for determination of electric-dipole matrix elements (or 
lifetimes) of principal transitions from these coefficients. 

The method is illustrated for atomic Cs. In recent 
years the most accurate lifetimes for alkali-metal atoms 
were derived using photoassociation spectroscopy of ul- 
tracold atomic samples reviewed in Refs. [? ? ]. Un- 
fortunately, this approach was not as successful for Cs 
because of peculiarities of molecular potentials of the Cs 
dimer [? ]. We consider an alternative method. Based 
on the van der Waals coefficient C§ deduced from ultra- 
cold collision data, the calculations reported here lead to 
the matrix elements for Cs with an uncertainty similar 
to the best direct lifetime measurements. We note that 
the determined high-precision dipole matrix element for 
6S1/2 ~ 6P1/2 transition in Cs [? ] is crucial for an in- 
terpretation of the most accurate measurement [? ] of 
atomic parity nonconservation (PNC). Atomic PNC [? 
? ] provides powerful constraints on possible extensions 
(e.g. extra Z-bosons) to the standard model of elemen- 
tary particles. 

For Cs, ? ] deduced C 6 = 6890(35)! using results 
of high-resolution Feshbach spectroscopy [? ] . Recently 
more Feshbach resonances due to higher angular mo- 
menta were identified by Chu and co-workers and this 
value has been constrained further to C 6 = 6859(25) [? 
]■ 

The dispersion coefficient Cq characterizes a second- 



1 Unless specified otherwise, atomic units h 
used throughout the paper. 
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order dipole response of atoms to molecular fields at large 
internuclear separations. It can be expressed as a quadra- 
ture of dynamic polarizability [? ] 
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Here \v) is an atomic ground state (6S1/2 for Cs), AEi 
are energies of intermediate states |i) taken with respect 
to the ground state, and D are electric-dipole matrix ele- 
ments. For Cs almost 85% of Cq is accumulated from 
intermediate states 6P1/2 and 6P 3 / 2 - We exploit this 
strong dependence to deduce matrix elements of princi- 
pal transitions by calculating residual contributions using 
ab initio methods. We separate the contribution of 6Pj 
states to dynamic polarizability and write 



a iw 



a p (iu) + 



a r (iu) 



where a r combines contributions of other intermediate 
states. These polarizabilities are shown in Fig. [[]. Intro- 
ducing reduced matrix elements Dj = (6Si/ 2 \\D\\6Pj) 
and a ratio R = (D 3 / 2 / 'D 1 / 2 ) 2 , 




(3) 

Energies AEj of 6P7 states arc known experimentally 
with a high accuracy and the ratio of matrix elements 
R = 1.9809(9) has been measured with a 0.05% precision 
by ? ]. This accuracy is much better than the errors 
in individual matrix elements due to a cancellation of 
systematic uncertainties. 
Since 



C fi = - 



du 



a p (iu) 2 + 2a p (iu) a r (iu) + a r (iu) 2 
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FIG. 1: Comparison of contributions to dynamic dipole po- 
larizability from the principal transitions a p (iiS) (dashed line) 
and residual intermediate states a r (iuj) ( long-dashed line). 
Total polarizability a (ito) is represented by a solid line. 
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FIG. 2: Comparison of deduced (filled circles) values for 
matrix elements {6^1 / 2 1 |-D| |6S , i/ 2 ) and {6P3/2WDW6S1/2) with 
the results of direct lifetime measurements [? (filled 
squares). The values are normalized to the present results 
and the deviation in % is shown. 



the van der Waals coefficient can be parameterized in 
terms of the matrix element -D1/2 as 



where 
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Solving the quadratic equation (|^) we obtain 
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and the problem is reduced to an accurate determination 
of parameters £. The calculation of these quantities and 
uncertainty estimates are discussed below. We find £ p = 
14.0787, £ x = 46.05(92), and £ r = 138.0(2.8) (the errors 
in the ratio R will be treated separately) and obtain 

<6P 1/a ||D||6S 1/a ) = 4.5028(60), (10) 
(6P 3/2 ||D||6S 1/2 ) = 6.3373(84). (11) 

These matrix elements are consistent with the results of 
direct lifetime measurements [? ] (6Pi/2||-D||65i/2) = 
4.4890(65) and (6P 3/2 ||.D|| 661/2) = 6.3265(77) and have 
a comparable accuracy (see Fig. ||). 

Details of calculation. - The parameter £ p can be 
evaluated with a high precision using experimental en- 
ergies and the ratio of matrix elements R from Ref. [? 
]. We obtain £ p = 14.0787. The uncertainty induced in 



this quantity by experimental error in the ratio R will be 
addressed separately. 

To determine parameters £ x and £ r we have to com- 
pute the residual dynamic polarizability a r {iw). Here we 
follow formalism laid out in Ref. [? ] and augment it with 
a refined error analysis. The intermediate states can be 
separated into two classes - valence states, both bound 
and continuum (a' v ) , and core-excited states (a c ) 

a r (iu>) = a v (iu>) + a c (iu>) + a cv (iu) . 

Here term a cv contains a small core- valence coupling cor- 
rection addressed below. 

In summation over valence states a' v we use a com- 
bination of relativistic linearized coupled-cluster singles- 
doubles (CCSD) method and Dirac-Hartree-Fock (DHF) 
approximation. In particular, the contribution from 7Pj 
and 8Pj states is computed with the matrix elements ob- 
tained in [? ] and experimental energies [? ] . The rest of 
the valence states is incorporated using B-spline quasis- 
pectrum [? ] generated in the "frozen" core DHF approx- 
imation. In the summation over core-excited states we 
employ relativistic random-phase approximation (RPA), 
described in Refs. [? ? ] with an obvious extension 
for frequency-dependence. In the RPA we allow for ex- 
citations of core electrons to all possible valence states, 
including the occupied 6S state. To account for a subse- 
quent violation of the Pauli exclusion principle, a counter 
core-valence coupling term a cv was introduced; we cal- 
culate it in the DHF approximation. 

Now we proceed to the calculation of parameter £ x . We 
notice that in Eq. (0) a smooth and broad curve a r (iu>) 
is integrated with a narrow Lorentzian-shaped function 
a p (iuj) (see Fig. |l|). To illustrate the main source of 
uncertainty we approximate 
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duj = a r (0) (1 + R) , (12) 



1. e. 



the uncertainty in £ x is governed by static residual 
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TABLE I: Breakdown of various ab initio contributions to 
static dipole polarizability from the intermediate states be- 
yond 6Pj. 





States 


Value 


Method 


a v (0) 


7Pj 


1.50 


CCSD" 




8Pj 


0.18 


CCSD" 




9 • • ■ ooPj 


0.33 


DHF 




Total 


2.01 




a cv (0) 




-0.47 


DHF 


Oc(0) 




15.81 


RRPA 


Q r (0) 




17.35 





a ab initio matrix elements and experimental energies. 



polarizability a r (0) . We present a breakdown of various 
contributions to this polarizability in Table || 

The error bars of the derived matrix elements depend 
sensitively on the uncertainty of £ x . To estimate this 
uncertainty, it is instructive to discuss characteristic ac- 
curacies of various relativistic many-body methods. As a 
test case we consider contribution of principal transitions 
to the static dipole polarizability 

We calculate this quantity in (i) Dirac-Hartree-Fock 
(DHF) approximation, (ii) linearized coupled-cluster 
singles-doubles (CCSD) method [? ], and (iii) using 
experimental energies [? ] and matrix elements [? ]. 
The results are presented in Table 0. It is clear that 
the Hartree-Fock approximation has a 70% error and the 
CCSD values have an accuracy in the order of 1%. Based 
on this example we assign a 1% uncertainty to CCSD 
contributions from 7Pj and 8Pj intermediate states and 
a 70% uncertainty to the corrections from the rest of 
the valence state and 70% to core- valence coupling term. 
We further replace the RPA value for the static polar- 
izability of the core with accurate semicmpirical value 
15.644(5) [? ] and obtain < c (0) = 17.18(40). Almost 
entire error comes from the uncertainty of the Hartree- 
Fock method. It is worth noting that the overall accu- 
racy of the a v (0) and the derived matrix elements can 
be further improved using, for example, linear-response 
coupled-cluster method [? ]. Based on Eq. ( |i~2| ) we im- 
prove the accuracy of the parameter £ x by rescaling the 
ab initio value of the integral, Eq. (^) , with the semiem- 
pirical a^.° (0). The final value for the quantity £ x is 
46.0(1.0) . 

The overwhelming contribution to parameter £ r comes 
from the dynamic core polarizability a c (iw). In sum- 
mation over core-excited states we employ relativis- 
tic random- phase approximation (RPA), described in 
Refs. [? ? ] with an obvious extension for frequency- 
dependence. Because of the equivalence of length- and 
velocity-forms of dipole matrix elements in RPA, the 
calculated dynamic polarizability of the core satisfies 
an important Thomas-Reiche-Kuhn (TRK) sum rule 



TABLE II: Characteristic accuracy of ab initio methods in 
calculations of contribution of principal transitions to static 
polarizability of Cs ground state. Accuracy of DHF and 
CCSD methods is defined with respect to the experimental 
value (third row). 



Method 


Value 


Accuracy 


DHF 


644 


70% 


CCSD" 


378.6 


0.7% 


Expt. 


381.2 


0.2% 6 



a ab initio matrix elements and experimental energies. 

b Based on experimental error bars in matrix elements from Ref. [? 

] 




FIG. 3: Sample Brueckner-Goldstone diagram illustrating 
coupling of a doubly excited core state (double horizontal 
line) to the many-body ground state of a closed-shell core 
through a one-body dipole matrix element (wavy line). Such 
contribution is possible only in a combination with correla- 
tion part of the wavefunction of the ground state (first order 
correction is shown, with Coulomb interaction represented by 
dashed line.) 



lim^^oo a c (iu>) = N c /u 2 , N c being the number of core 
electrons (N c = 54 for Cs.) This property is illustrated 
in Fig. 1 of Ref. [? ]. The static (uj = 0) core polar- 
izability obtained with RPA is a c (0) = 15.81. This re- 
sult is in a 1% agrement with a value of 15.644(5) de- 
duced from semiempirical fitting of Rydberg spectrum of 
Cs [? ]. The RPA accounts only for a correlated con- 
tribution of particle-hole excitations. However, the cor- 
rection from multiple core excitations is expected to be 
strongly suppressed because they can couple only to the 
first-order many-body correction to the core wavefunc- 
tion (see Fig. ||). In addition, the effect of multiple core 
excitations is reduced by larger energy denominators in 
the expression for polarizability. 

Overall the RPA approximation results in a dynamic 
core polarizability differing from the high-precision value 
by 1% at to = and at large frequencies recovering the 
exact TRK limit [? ]. The parameter £ r calculated in- 
cluding only core polarizability, £ r = 137.2, is slightly 
increased by inclusion of valence states and the core- 
valence coupling term a cv to 138.0. Based on the above 
discussion we assign a 2% uncertainty to £ r and arrive 
at = 138.0(2.8). It should be noted that more conser- 
vative error bar of 5% for £ r does not significantly affect 
the resulting precision of derived matrix elements. 

Employing C§ — 6859(25) [? ] and the calculated pa- 
rameters £ we determine the matrix element Di/ 2 , and, 
using the ratio R [? ], the matrix element D-i/i- The 
obtained values are given in Eq. dl(i|). The uncertainty 
in matrix element calculated from Eq. (j9j) can be param- 
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eterized as 



SD 
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D 1/2 
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where A Ce = 7.5 x 1CT 2 , A = 3.0 x 1CT 5 , Ar = 1.4 x 
10~ 3 , and An = 9.4 x 10~ 2 . Combining estimated un- 
certainties we obtain an error bound of 0.13% for matrix 
elements. This accuracy is similar to that of the best 
direct lifetime measurements by ? ]. The reader is re- 
ferred to Ref. [? ] for an extensive comparison with 
other measurements and theoretical predictions. Our re- 
sult is most sensitive to the errors in the van der Waals 
coefficient and parameter £ x . Provided that Cg is known 
exactly, the uncertainty in the matrix elements can be re- 
duced to 0.07% with the current technique for estimation 
of £ x . It is worth noting that the precision of calculation 
for £ x can be substantially improved, for example, using 
linear-response coupled-cluster method [? ]. 

Conclusion. — We exploited a strong dependence of the 
van der Waals coefficient Cq on matrix elements of prin- 
cipal transitions. We deduced these matrix elements by 
calculating small residual contributions using ab initio 
methods. The proposed method was applied to Cs atom 
and the derived matrix elements are consistent with the 



best direct lifetime measurements [? ] and have a similar 
uncertainty. 

In anticipation of high-precision measurements of 
static dipole polarizability a(0) of the ground state of Cs 
with atom interferometry [? ], we note that 96% of the 
polarizability is due to the contribution of the principal 
transitions a p (0) [? ]. Subtracting the residual contribu- 
tion a£' e '(0) = 17.18(40) from the measured a(0) one can 
also determine matrix elements of principal transitions to 
an accuracy of 0.05% from Eq. (g). 

We employed the 0.36%-accurate value of the van der 
Waals coefficient Cq deduced by ? ] from high-resolution 
Feshbach spectroscopy of ultracold Cs atoms [? ]. Ac- 
cording to Ref. [? ] the precision of this dispersion co- 
efficient can be potentially improved to 0.03%. If such a 
precision is achieved, the method proposed here, aug- 
mented with more accurate ab initio calculations and 
better measurements of the ratio R, can lead to deter- 
mination of matrix elements of principal transitions with 
an unprecedented 0.01% accuracy. 
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